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The application of methods of time-dependent density functional theory (TDDFT) to systems of 
qubits provided the interesting possibility of simulating an assigned Hamiltonian evolution by means 
of an auxiliary Hamiltonian having different two-qubit interactions and hence a possibly simpler wave 
function evolution. In this note we extend these methods to some instances of Lindblad evolution 
of a spin chain. 
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I. INTRODUCTION 

The application of methods of time-dependent density 
functional theory to systems of qubits has been proposed 
by Aspuru-Guzik and Tempel in reference In particu- 
lar, these authors have given, for a family of Hamiltonian 
evolutions of qubit systems, a constructive proof of the 
van Leeuwen theorem (VL) Q. 

We review their work on the Hamiltonian evolution of a 
spin chain in section II: we explore the possibility offered 
by van Leeuwen's mapping from densities to potentials 
of simulating, by means of a spin chain having spatially 
homogeneous hopping parameters, the magnetization of 
an engineered spin chain Q. 

Section III explores the possibility of extending these 
methods to some instances of dissipative evolution of a 
spin chain according to the Lindblad equation. We focus 
our attention on a class of models studied in all details 
by Znidaric 0(3 and examplified in figure Q] 
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FIG. 1. The end points of a spin chain of length s interact 
with two reservoirs Bl, Br. Dephasing operators act inde- 
pendently on each site. 

The case in which only dephasing is present is of some 
interest on its own right: in subsection III. A we numeri- 
cally explore the possibility (suggested in for electronic 
systems) of simulating, in this case, the magnetization of 
an open quantum system with that of a driven closed 
system; then we discuss the complementary problem of 



compensating for the effects of dephasing on a spin chain. 
Section III.B is devoted to the study of the dissipation 
current Q @ present when we allow exchange of magne- 
tization with the two reservoirs Br. 

In section IV we discuss the difficulties, already evidenced 
in @, that appear when the local kinetic energy is not 
everywhere nonzero, and present suggestions for further 
work. 



II. THE MODELS: HAMILTONIAN 
EVOLUTION 

We consider a chain of s spin 1/2 systems 
<x(l),<x(2),...,er(s), with a(x) = (ai(x), a 2 (x), <r 3 (x)), 
evolving under the action of a Hamiltonian of the form: 



H(t) = J2 J{x){a x {x)a 1 {x + 1) + a 2 (x)a 2 (x + 1)) + 

x=l 

s— 1 s 

+ ^2 K(x)a 3 (x)a 3 (x + 1) + ^ h(x, t)a 3 (x) = 

X — l X—l 

s 

= H (J,K) + J2h(x,t)a 3 (x). (1) 

x=l 

As our notations indicate, the two-qubit coupling con- 
stants J(x) and K(x) are supposed to be time indepen- 
dent, while the field h(x, t) is allowed to depend on time. 
The current operator field j is defined by: 

j(0) =j(s) = 

j(x) = 2J(x)(a 1 (x)a 1 (x + 1) - a 2 (x)a 2 (x + 1)). (2) 

Its interest comes from the following commutator identity 
(from which the continuity equation easily follows): 

- i[a 3 (x),H(t)] + (j(a;) - i(x - 1)) = 0. (3) 

It is useful, in what follows, to consider also the identity: 

= 4iJ{x)(<7+(x)a-(x + 1) - a_(x)a+{x + 1)), (4) 



* defalco@dsi.unimi.it, tamascelli@dsi.unimi.it 



where a±(x) = {<ti{x) ± i<r 2 (x)) /2. 

Using the above identities it is immediate to prove the 
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following further identity, of obvious relevance to the 
Heisenbcrg evolution of the current field: 



(5) 



-i\i(x),H(t)} = 

8J(x) (J(x) (a 3 (x) - a 3 (x + 1)) + 
+ J(x - l)a 3 (x)r(x -1,x+ 1)+ 

- J(x + i)a 3 (x + 1)t(x, x + 2)+ 
+ K(x + l)a 3 (x + 2)t(x, x + 1)+ 

- K(x - l)o- 3 (a; - l)r(x, x + 1)+ 
+(h(x + 1, t) - h(x, t))r(x, x + 1)) . 

Here and in what follows we set 

r(x,y) = a + (x)(T-(y) + a-(x)a+{y). 



Identities ((3]) and j5]) have the obvious consequence that, 
if the state | ip(t) ) satisfies the Schrodinger equation 



i-\m)=H(t)\m), 

then its mean magnetization field 

m 3 (x,t) = (ij(t)\o- 3 (x)\ip(t)) 
satisfies the continuity equation 
dm 3 (x, t) 



(0) 



dt 



+ (j(x,t)-j(x-l,t)) = 



(7) 



and its mean current field 

satisfies the equation of motion 
d 



dt 



j(x,t) = \[i(x),H (J,K)]\ m )- 



(8) 



+ 8J(x)(h(x + l,t)- h(x,t))( ip(t) \t(x,x + 1)| tfj(t) ). 

In reference [TJ Tempel and Aspuru-Guzik take the fol- 
lowing point of view: suppose we are interested in attain- 
ing, in a neighborhood of t = 0, a preassigned magnetiza- 
tion field m 3 (x,t); because of ([7]) and ((2]) this target will 
be reached if we are able to attain a current field j(x, t) 
given by 



j(x,t) = -^^2,m 3 {x,t); 

y<x 



(9) 



finding a Hamiltonian Hq(J,K) + ^ s x= ih(x,t)a 3 {x) 
(where J and K are supposed to be assigned, and h is 
supposed to be unknown) that determines such a cur- 
rent field j(x,t) is, in turn, equivalent to solving the 
differential- algebraic problem 

jTj{x,t) = -i(^(t) \\i(x),H (J,K)}\ i,(t)) + 
8J(x)(h(x + l,t) - h(x,t))(ip(t) \t(x,x + 1)| -0(t) ) 

[ ii\m) = (MJ,K)+Y, S x =iHx,t)a 3 (x))\i>(t)) 

(10) 



in the unknown force field h(x + 1, t) — h(x, t) and in the 
unknown wave function \ ip(t) ). Naturally, | ip(t) ) is rc_ 
quired to have the correct initial magnetization m 3 {x, 0) 
and the correct initial current j(x, 0). For this formu- 
lation of the van Leeuwen theorem [2| as a non linear 
Schrodinger equation in which the potential is itself a 
functional of the wave function, we refer the reader to 
references !] and [U. 

In the applications, Hq(J,K) typically describes the 
two-qubit interaction one is physically able to achieve. 
Suppose however that the task one has to undertake 
is the transfer of magnetization along the chain [ll| 
[JJJ realizable through a different two-qubit Hamilto- 
nian Hq(J,K). The question is: can a force field 
h(x+ l,t) — h(x, t) compensate for the "wrong" choice of 
two-body couplings? 

As an example, consider as a target system a spin chain 
evolving under the engineered XY- Hamiltonian: 



^ 1 
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~K\J x(s — x) 



t(x, X + 1). 



(11) 



The Hamiltonian (fTTj) . first discussed in [13[ and then 
rechristened engineered XY-chain in 0], is quite interest- 
ing for the quantum computing community: it induces a 
periodic behavior of the system with period T = 2s; in 
particular it realizes the perfect transfer of an excitation 
(spin up) initially located at the first site x = 1 to the 
final site x — s in a time s fl2l ]. 

Figure [2] summarizes some experience we have gained in 
the numerical integration of (|10p for s = 6, starting from 
the initial condition 



(12) 



where we have indicated by 

\x} = \ cr 3 (x) = +1, a 3 (y) = -1 for y ^ x) 

the simultaneous eigenstate of the operators a 3 in which 
only the spin at position x is "up " . 
As figure [2] shows, without undertaking the difficult task 
of implementing the engineered couplings, varying as 
x (s — a;), it is possible to obtain the same magnetiza- 
tion and current fields using spatially homogeneous two- 
body couplings and suitably chosen control fields h c (x,t), 
e.g. by letting, as we did, the initial condition |^ ) evolve 
under 
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+ 1) + yh c (x,t)a 3 (x). (13) 



x=l 



III. OPEN SPIN CHAINS 

We discuss the possibility of extending the above dis- 
cussion to the case in which the initial state of the system, 
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described by a density matrix poj evolves according to a 
Lindblad equation of the form: 



dp(t) 
dt 



i[p{t),H{t)]+V{p(i)) 



(14) 



where H(t) is of the form (JTJ). 

The dissipator V that we consider takes into account 
both the dcphasing induced by the interaction of each 
spin with some external degrees of freedom (for example 
phonons) and the coupling with unequal one-spin baths 
[lj] Bl and Br at both ends of the chain; it has the form: 



V(p) = V deph (p)+V bath (p); 



(15) 



The dephasing superoperator is given, in Lindblad form 
HI, by 



V deph (p) = £ {[L(x)p,L(x)l] + [L(x),pL(x?]) 



x=l 



(16) 



where the generator associated with each site x is 



L ( x ) = \l ~^ a 3( x )- 



V , instead, involves four generators: 



L x = Ve(l-MV+(l), L 2 = vMl + M)a_(l), (17) 



L 3 = \/e(l+p)o- + (s), L 4 = y/ e(l - /x)cr_(s), 

where the coupling parameter e and the asymmetry pa- 
rameter p satisfy the conditions e > and — 1 < p < 1. 
The asymmetry parameter /j models a possible difference 
of the chemical potentials of Bl and Br. We refer the 
reader to the second section of [HI for a discussion of 
the range of validity of this approach and of the possible 
degree of control on the baths and on the bath-system 
interaction. 

For the sake of clarity, we discuss the effect of the two 
sources of dissipation separately. 

If there is no interaction with the baths (e = 0), it is, for 
1 < x < s: 

Tr(a 3 (x)V(p))=0. (18) 
The continuity equation ([3]) retains, then, the form 
dm 3 (x, t) 



dt 



+ (j(x,t)-j(x-l,t)) = 



where it is, now, 7773(2;, i) = Tr (o~ 3 (x)p(t)) and 
j(x,t)=Tr(](x)p(t)) : 

Simple algebraic manipulations show that, for 1 < x < s, 
it is: 



Tr(j(x)V(p))=-4rfrr(j(x)p). 



(19) 



The evolution of the current field j(x,t) is modified by 
the presence of dcphasing, namely: 



dj{x,t) 
dt 



= -i Tr (\j(x), H(t)]p(t)) - 4 V j(x, t). (20) 



The commutator on the right hand side is given by (f5|). 

The evolution in the presence of the interaction with the 
baths Bl and Br (77 = 0, e > 0), on the other side, de- 
termines a nontrivial modification of the continuity equa- 
tion. For 1 < k < 4, we set 



VL k (p{t)) = [L k p,L\ 
so that it is 



L k ,pLi 



V bath (p) = '£'DL k (p). 



(21) 



k=l 



A little algebra shows now that the effect of the left bath 
Bl on the magnetization is: 

Tr (a 3 (x) (VL 1 (p)+VL 2 (p))) = 
= -6 1>x 4e(p + Tr(a 3 (x)p)) 

whereas the effect of the right bath Br is: 

2V(<7 3 (aO (M3W+M4W)) = 
= -S StX Ae(-p + Tr(a 3 (x)p)) . 

The continuity equation reads now: 

dm 3 (x, t) 



dt 



+ (j{x,t)-j{x-l,t)) 



(22) 



= -4e (Sx !X (m 3 (x, t) + p) + S S!X (7773(2;, t) - p)) . 

The rhs of (|22| plays the role of a dissipation current [3] ■ 
As to the time evolution of the current field it is: 



dj{x,t) 
dt 



■■-i Tr(\i(x),H(t)]p(t))- 
2e (S l x + 6 s _ l x )j(x,t). 



(23) 



Summarizing, the kinematical and dynamical equations 
in the presence of both dcphasing and coupling with the 
baths (77 > and e > 0) are: 



dm 3 (x, t) 
Jt 



(j(x,t)-j(x-l,t)) 



(24) 



= -4e (Sx^ (7773(2;, t)+ p)+ S S:X (7773(2;, t) - p)) 

d ^ = -i Tr(\j(x),H(t)]p(t)) + 
- (477 + (S hx + 5 s - hx ) 2e) j(x, t). 

We conclude this section with a few numerical examples. 
In all the examples that follow, we take an initial condi- 
tion of the form po = \ ip$ ) ( ip$ \ , with | -ipg ) given, for 
different values of s, by ((T2^) . 
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A. Pure dephasing 

In the case of pure dephasing, the observable 

l + cr 3 {x) 



N 3 



E 

x = l 



namely the number of particles (spins "up"), is a con- 
stant of motion. 

It has been shown in reference @ that it is possible to 
reproduce current and density of an electronic system 
subject to a number-conserving dissipation by letting 
the same quantum system evolve unitarily under a time- 
dependent Hamiltonian with suitably chosen control field 
h c (x,t). 

We wish to show that the same considerations apply to 
open spin chains as well. 

Suppose that the interactions in the chain are described 
by the Ay-Hamiltonian: 

1 s_1 

H " = -o £ M x )°-( x + !) + °+(* + !)"-(*)) (25) 

x=l 

that can be obtained from the Hamiltonian family de- 
scribed by (TI|) by setting: J(x) = —1/4, K(x) — 
h(x,t) = 0. 

The evolution of the chain under the sole effect of de- 
phasing from the initial state po is determined by the 
Lindblad equation dHJ) with V = V deph defined as in 
(|16[). We indicate by j deph (x,t) and m 3 eph (x , t) the en- 
suing current and magnetization fields. 
We wish to determine an external field h c (x, t) such that 
the evolution p c (t) of the initial state po under the Hamil- 
tonian: 



H c {t) =H + J2h c (x,t)<j 3 (x) 



(26) 



determines the same magnetization and current fields as 
in the dissipating system, that is: 

m c 3 (x, t) = Tr (a 3 (x)p c (t)) = mf ph (x, t) 
f(x,t)=Tr(i(x)p c (t))=j d ^ h (x,t). (27) 

It is sufficient to solve, for the state p c {t), the Liouville 
problem 



P c (0) = Po 
d P c (t) 



dt 



i\p c {t),H c {t)} 



(28) 



where the unknown field h c (x,t) is itself a functional of 
p c through the algebraic equation: 



dj deph (x,t) 



dt 



= -i Tr([j(x),H }p c (t)) + 



(29) 



In the numerical example that we discuss below, 
dj deph (x,t)/dt is known, having been determined by nu- 
merical integration of the Lindblad equation for pure de- 
phasing under the initial condition | V'o ) ( V'o I with s = 20 
and r\ = 0.01. 

Figure [3] shows that, for this example, the differential- 
algebraic problem posed by (|28[) and (|29p does have a 
solution that very well attains the current of the dissi- 
pating system (j c = j deph ) and, because of the continuity 



equation also the magnetization (m 3 



deph 



The Kohn-Sham potential |6j shown in frame 3. (a) is 
able to reduce the amplitude of the oscillations visible, 
for Hamiltonian evolution under Hq, in frame 3.(b); this 
leads to the magnetization profile of frame 3.(c) in which 
the controlled magnetization is reported: it is numer- 
ically indistinguishable from the target magnetization 
(not reported). Because of the left- right symmetry of 
our problem we have found it sufficient to let x go only 
from 1 to s/2. 

As a further example, we pose the following, complemen- 
tary, problem: we are given an open quantum system 
affected only by dephasing (e = 0, T) > 0) and evolving 
under a master equation of the form (|14|) . Is it possible 
to determine an external control field h c that compen- 
sates for the effects of dissipation, at least as far as the 
magnetization m 3 and the current j fields are concerned? 
For the sake of definitcness, we consider a spin chain with 
Ay-interaction Ho, as in (125[) . The system is affected by 
noise: the dissipator V deph (see equation (Tj"6|) ) intervenes 
in the evolution of the system as in (fT4^) . 
We are looking for a choice of a control field h c such that 
the current and magnetization fields evolve in time as if 
the dissipator were not acting at all; namely the target 
is: 

m 3 (t,x) = Tr (cr 3 (x)p(t)) 
](x,t)=Tr(j(x)p(t)) 
where ~p(t) is the solution of: 



f t =i[H,p(t)] 



Po- 



+ 8J(x) (h c (x + 1, t) - h c (x, t)) Tr (t(x, x + l)p c (t)) 



with the initial condition p(0) 

Figure 3] refers to a system of s = 3 spins. It shows 
that our control scheme works only locally; in fact, as 
soon as one of the local kinetic energy terms T{x) = 
2J(x)t(x,x + 1) hits the value zero, the corresponding 
control field h c develops a singularity, as shown in figure 
Ha). 

It is not surprising that the control problem posed above 
admits only a solution that is local in time. This as pec t 
of the problem has already been discussed in Q and [l6l | . 



B. Dissipation current 

In order to discuss the role of equation (|22p . we ana- 
lyze here, for the sake of definitcness, the same problem 
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posed at the end of the previous subsection (is there an 
external field that "compensates" for dissipation?) in the 
presence of interaction with the baths Bl, Br. 
The novel feature that appears in the case e > is the 
fact that the number operator iV 3 ceases to be a con- 
stant of motion: "charges" in the form of "spins up" can 
be pumped into or drained from the chain by the baths. 
Because of this fact, one must carefully distinguish be- 
tween magnetization and current. In the discussion that 
follows, we set 

s 

2 = 1 

with Ho given by equation (f2"5j) . and consider the dissi- 
pator flU) 



fe=i 

for the definition of the current operator j(x) we refer 
to equation ^ and, as usual, we take in our numerical 
example p = \ ijj ) ( ipo I with \tpo) as m (EJ) • 
Figure [5] answers, in the particular case s = 3, the follow- 
ing question: is there a field h c (x 1 t) such that, under the 
action of a Hamiltonian of the form H c (t), and of the dis- 
sipator T> bath (p), the initial condition p(0) evolves into a 
state p c (t) for which the expectation value Tr (p c (t)j(x)) 
coincides, for 1 < x < s — 1 and at least in a neighborhood 
of t = 0, with 

— def 

j(x,t) = (ipo\cxp(+itH )i(x)exp(~itH )\'il}o)? 

The behaviour of the control field h c (x,t) satisfying the 
condition 



Tr(p c (t)j(x)) = j(x,t) 



posed above shows (see figure [Ha)) that the existence of 
the solution is again local, because of the vanishing (fig- 
ure OJb)) of the expected local kinetic energy T(l) in the 
controlled evolution. 

Figures [5jc) and[Hd) show that, as long as it is defined, 
the field h c (x,t) does enforce the target current: this is 
hardly surprising since it has been found by imposing 
precisely this condition. 

In order to compare, in the figures, the controlled evo- 
lution with the Lindblad evolution without control, we 
define p actua( (£) as the evolution of the initial condition 
Po determined solely by the Hamiltonian Hq and the dis- 
sipator V bath . 

Setting 

m 3 (x,t) = ( ipo | exp(+itH )<j 3 (x) exp(-itH )\ ip ), 

figure [S] answers, in the same situation as above, the fol- 
lowing different question: for which values of x does the 



state p c (t) determined by imposing the condition (|30p 
satisfy the additional condition 



Tr (p c {t)<j 3 (x)) =fn 3 (x,t)? 



(31) 



In order to answer this question, we must use equation 
(1221). For 1 < x < s it is: 



dm 3 (x, t) 
Jt 



(f(x,t)-f(x-l,t)) = 
dm 3 (x, t) 



and 



(j{x,t)-j(x-l,t)) 



m 3 (x, 0) = 77i3(x, 0). 



dt 



Therefore, the controlled magnetization m 3 (x,t) at sites 
2, 3, . . . , s— 1 is, as long as it exists, identical to the target 
one (figure [U(a)). 

As to the sites coupled to the reservoirs, because, again, 
of equation (|22p , the following differential equation holds 
for m 3 (l, t) 



dm c 3 (l,t) 
Jt 



-f(l,t)~4ep-4em c 3 (l,t) 



-j{l,t)-Aefi-4em c 3 (l,t) 
dm 3 (l, t) 



dt 



Aep - 4em§(l,i). 



Under the initial condition m§(l,0) = m 3 (l, 0), this im- 
plies that 

m{j(M) = m 3 (l,t)-/i(l-e- 4et ) + 

-4ee- 4e * f m 3 (l,T)e 4£t dr. (32) 
Jo 



(30) Similarly, a site s, it is 



m c 3 (s,0) = m 3 (s,0) 
dm 3 {s,t) drn 3 (s,t) 



dt 

and, therefore 



dt 



+ 4e/i - Aem c 3 (s,t) 



m c 3 (s, t) = m 3 (s, t) + ii (1 - e~ 4£t ) + 

-4ee- 4e */ m 3 {s,T)e 4et dT. (33) 
Jo 

Therefore, the controlled magnetizations at sites 1 and s 
are, as long as they exist, different from the target ones 
(see the initial part of figure E^b)). 

By explicit computation it is easy to see that, in the 
simple example s = 3 at hand, the target magnetizations 
are given by 



m 3 



(2,t) = -^cos (y2tj 



6 



and 



m 3 (l,t) = m 3 (3,t) 



6 



V2t 



With the above explicit expressions of 777,3 and m§(3, t) it 
is immediate to compute the right hand side of equations 
{32) and {33]), namely 



a(t) = 7773(1, t)-n(l- e- 4ct ) - 4ee~ 4£t / m 3 (l, r)e 4et dr 
6(t) = m 3 (s, i) + /1 (1 - e' 4et ) - 4ee- 4et / 7tT 3 (s, r)e 4£t dr. 

These are the two quantities represented by the thin 
dashed and dotted lines in figure [6jb), extending well 
beyond the interval < t < 3.55 of existence of the solu- 
tion of the VL problem. 

The fact that there is a value of t beyond which a(t) and 
b(t) can go out of interval [—1, 1] of physically acceptable 
values of magnetization, convinces us of the fact that the 
determination of a control field compensating for the dis- 
sipation induced by the interaction with the baths Bl and 
Br admits a solution that is at most local in t. The local 
character of the solution is therefore determined neither 
by the choice of the initial condition, nor by the Hamil- 
tonian and dissipator, nor by the numerical algorithm of 
solution of the differential algebraic problem. 



IV. CONCLUSIONS AND OUTLOOK 

We have probed the interesting idea, advanced in ref- 
erence [l| , of quantum computing, by means of spin 1/2 
arrays, without wavefunction. 

In our exposition we have stressed the dynamic con- 
trol perspective advanced in [l[: the components of 
the wave function or, in our generalization to the case 
of open quantum systems, of the density matrix, are 
the latent dynamical variables, that in their evolution 
(parametrized by J, K, e, \i, 77, | ip® )) determine the re- 
sponse, in terms of output magnetization and/or current, 



to the time-dependent input variables h c (x,t). 

In our choice of numerical examples, we have privileged 
themes motivated by our previous experience with Feyn- 
man's model of a quantum computer 
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For instance, the numerical example of figure [5] can be 
rephrased in the following way: can one use the free- 
dom of reparametrizing the system (J, K, . . . , \ tpo ) ~ * 
J', K ', . . . , I i/'q )) without altering the input- output re- 
sponse in such a way as to simulate, on a transitionally 
invariant structure, the delicate tailoring of hopping pa- 
rameters required for perfect transfer of the cursor along 
the program line? 

The example of figure [4] can be similarly reformulated as: 
can one compensate by external fields for, in the words 
of Feynman jl7| , imperfections and free energy losses in 
the motion of the clocking cursor! 

In answering the above questions, the success of TDDFT 
has been only partial because the solutions exist only lo- 
cally in time: the point is that in equations (jTUJ) and (f2"5)) 
the unknown gradient h c (x + 1, t) — h c (x, t) appears with 
the local mean kinetic energy as a coefficient (a handle). 
This difficulty would turn out to be particularly severe if 
one tried simulate a sharp initial condition in the Hamil- 
tonian context of figure [2] the handles are in this case 
initially zero, to start with. 

In the context of figure 2] (dephasing), vanishing of the 
handles is unavoidable, because they are precisely the 
nondiagonal matrix elements of the density matrix that 
dephasing is supposed to damp. 

The example of figure [5] and figure [5] (interaction with 
baths) shows that it is indeed possible to compensate 
for dissipation even when the number of particles is no 
longer a constant of motion, at least for the spins that 
are not coupled with the baths. It is not surprising that 
the magnetizations of the extremal sites is not attainable 
by means of the control field: the field h c does commute 
with the number operator and therefore is not able to 
absorb injected charges or to create new ones. 

Further efforts are needed, we think, in the choice and/or 
in the smoothing of the initial condition and in the in- 
troduction of "convergence factors" [l[ that extend the 
existence interval of the solution with little errors in the 
overall propagation. 
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FIG. 2. s = 6; Frame (a): the magnetization field determined by the target Hamiltonian Frame (b): the magnetization 

field determined by the controlled Hamiltonian (|13[) . Frame (c): the control field h c (x,t), obtained by numerical integration of 
(|fO|l . as a function of x and t; interpolation in x is only for graphical convenience. Frame(d): the control field h c (2,t) (solid 
line) and the the mean value of r(l,2) in the controlled state (dashed line). In order to avoid singularities in h c (x,t) due to 
the zeroes of this mean value we have used Tikhonov regularization 




FIG. 3. s = 20, r) — 0.01, e — 0. Frame (a): the Kohn-Sham potential h c (x,t). Frame (b): the field 7713(2;, t) corresponding to 
purely unitary evolution of the initial condition po under the Hamiltonian (125[) . Frame (c): the magnetization field m^(x,t). 
Frame (d): 7713(10, t) (thin line) and 7713(10, t) (thick line). The field mf :ph (x,t) is numerically indistinguishable from m§(a;,t). 
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FIG. 4. s = 3, r\ — 0.01, e = 0. Frame (a): the control field h c (2,t) is represented (rescaled by a factor 0.1 for graphical 
convenience) as a function of time as a dot-dashed thick line. Expectations of the local kinetic energy T(l), under three 
different evolutions, are represented as functions of time: the dashed line corresponds to pure Hamiltonian; the solid thick line 
to Hamiltonian+dephasing; the solid thin line to Hamiltonian+dephasing+control field. Frame (b): the expectations of 03(1) 
for the same three systems. Same graphical conventions as in frame (a); the dashed line is not visible since the magnetizations 
of the target and controlled systems are, as long as both exist, numerically indistinguishable. 
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(c) (d) 



FIG. 5. s = 3, e = 0.1, fi = -1, 7/ = 0. Frame (a): the field h c (x,t). Frame (b): Tr (T(l)p(i)) (dotted line); Tr {T{l)p actua \t)) 
(dashed line); Tr (T(l)p c (t)) (solid thick line). Frame (c): j(l,t) (dashed line); f{l,t) (solid thick line); j actual (l,t) (solid thin 
line). Frame (d): j(2,t) (dashed line); f(2,t) (solid thick line); j ac ^ ual (2,t) (solid thin line). The current field j c (x,t) is, as 
long as it exists, numerically indistinguishable from the target field j(x,t). 
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(a) (b) 



FIG. 6. s = 3, e = 0.1, fx = -1, 77 = 0. Frame (a): Tr (0- 3 (2)p ac *" a! (t)) (solid thin line); m 3 (2,t) (dashed line); m§(2,t) (solid 
thick line). The magnetization m c (2,t) is, as long as it exists, numerically indistinguishable from the target magnetization 
7713(2;, t). Frame (b): a(t) and b(t) as functions of time (dashed and dotted thin lines, respectively); the corresponding fields 
7713(1, i) and 772.3(3, t) determined by numerical (local) solution of the VL problem are represented as thick lines; solid thin line: 
the target magnetization 7773(1,!;) = 7773(3, t). For the extremal sites of the chain, the controlled system does not attain, even 
locally, the target magnetization. 



